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We study a homogeneous Bose gas with purely repulsive forces. 
Using the Kac scaling of the binary potential we derive analytically 
the form of the thermodynamic functions of the gas for small but finite 
values of the scaling parameter in the low density regime. In this way 
we determine dominant corrections to the mean-field theory. It turns 
out that repulsive forces increase the pressure at fixed density and 
decrease the density at given chemical potential (the temperature is 
kept constant). They also flatten the Bose momentum distribution. 
However, the present analysis cannot be extended to the region where 
the mean-field theory predicts the appearence of condensate. 



1 Introduction 

Whereas the mean-field theory of an interacting Bose gas is now well under- 
stood, going beyond the mean-field description in a systematic way presents 
still an open problem (see jl] and references given therein). A natural tool 
for such a study is provided by the Kac scaling of the binary potential V{r). 
The scaled potential 
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becomes weak and long-range when the dimensionless parameter 7 approaches 
zero. However, the integrated potential energy 

a = J drV{r) = J drV^{r) < 00 (2) 

is 7-independent, and thus remains constant. The so-called van der Waals 
limit (i.e. the thermodynamic limit followed by 7 ^ 0) permits then to derive 
the mean-field effects corresponding to the potential V{r). They depend on 
the mean potential energy ap, where p is the number density of particles. 
For example, it has been shown that in the case of repulsive forces the van 
der Waals limit of the free energy density / yields the mean-field formula 

fmf = fo + (3) 

where /o refers to the reference system without interaction V{r). The formula 
(jni) holds both for classical and for quantum gases [21,|SI- In the theory of 
classical fluids one could go even further: the corrections to the zero order 
mean-field results have been derived, exhibiting the role of fluctuations for 
small but flnite values of 7 |lj . 

In the present paper we perform an analogous analysis for a homogeneous 
Bose gas, continuing our previous study p. Our object is thus to determine 
and to investigate the leading corrections to the mean-field theory predictions 
for equilibrium properties of a gas composed of identical bosons interacting 
via purely repulsive forces. The rigorous part of our analysis will be restricted 
the region of thermodynamic parameters in which quantum Mayer series 
(virial expansions) converge uniformly with respect to the small parameter 
7. In fact, under the assumption that both the potential and its Fourier 
transform are non-negative 

V{r) > 0, V{k) = j dr exp(-zk ■ v)V{r) > 0, (4) 

we can directly apply here the main results of fP. In fT] the Bose gas is repre- 
sented as a classical-like system of random polymers by using the Feynman- 
Kac path integral formulation of the Gibbs weight together with a cycle 
decomposition of the permutation group. Applying standard Mayer graph 
summation technique it has been shown therein that the density p and the 
chemical potential p of the Bose gas were related by a self-consistent relation 

p{p) = F{p - ap{p)) (5) 

The function F is defined by a convergent series of multiconnected graphs 
whenever the chemical potential is sufficiently negative, and the convergence 
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is uniform with respect to the scahng parameter 7. In this regime, and for 
7 <C; 1 the function F for a scaled potential (H)) takes the asymptotic form 

F^iu) = + 7^F(i)(z/) + 0(7^) (6) 

Here F/°gg(z/) represents the sum of zero order contributions from the tree 
diagrams. It turns out that 

F^rLi^) = Po(^) (7) 

where po denotes the perfect gas density. The limit 7^0 yields thus the 
self-consistent mean field equation 

Pm/(yu) = Po(/^ - aPm/(/i)) (8) 

Equation © implies that the dominant corrections to equation (jHJ are of the 
order 7^. The correcting term 

F^'\i^) = F,^rL{i^) + Fi2^iu) (9) 

results from the summation of the tree diagrams and from the summation of 
the ring diagrams. In Section 2 the form of the functions F^^^^lu) and F^llg{i>) 
is derived. The corresponding equation of state of the Bose gas is presented 
and analyzed in Section 3. It is remarkable that the dominant correction to 
the pressure has the same structure as that found by Hemmer for a classical 
fluid j3] . Under the additional assumption that the formula remains valid up 
to a neighborhood of the critical free density, we observe that the mean-field 
theory always underestimates the pressure in this region. We study then the 
one-particle reduced density operator displaying the nature of the momentum 
distribution (Section 4). Concluding remarks are presented in Section 5. 



2 Self-consistent equation beyond mean-field 

The function 

when used in the self-consistent equation (jH)) yields the term corresponding 
to the summation of the tree diagrams (section V.A in [Ij) . Here 

n^iq) =< expi-pU^) >, (11) 
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is the partition function of a single closed g— polymer, (3 = l/fc^T, T is the 
temperature, /c^ the Boltzmann constant, q the number of particles in the 
polymer and 

A = hJ^ (12) 
V m 

is the de Broglie thermal wave length. The asymptotic 7-expansion of the 
scaled self-energy f/^ of the polymer follows from the definition given in 

[/, = 7'^i^nO) + o(7^), 7-0 (13) 



,3^(^-1) 



SO that 

K-riq) = 1 - 7^m_^^v^(o) + 0(7^) (14) 
Inserting (fT^ into (fTUj) we get 

Ftleeil^) = Po{y) + l'Fi^l{u) + 0(7') (15) 

where 



V] 



7(1) 

tree\' / 2 



(16) 



We have adopted the notation ^ 



1 ^^/c 1 00 ^ 

/"V) = j.^M'^) = E ^exp(;3,.) (17) 



The dominant contribution from resummation of the ring diagrams 7 F, 



3f(1) 



ring \ 



has been also calculated in [T] (see eq. (70); it is derived here in the appendix 
for the sake of convenience, see eq. ()84|l ). It reads 

.(1) ^.^_if(2)^,^ f W{k)ff'\u) 



(27r)3 l + /5f(fc)/«(z/) 



Using the definition (fT7|) and denoting respectively by pg and by Pq the first- 
and the second-order derivative of the perfect gas density with respect to the 
chemical potential, we rewrite equations (fTH|l . (fTK|) in the form 

= ^Pi(W - ^Pi» (19) 



""^In ^ /j'^-' (z/) designates the 7-dependent function, including the vertex contribution 
K^{q). Here we have set 7 = 0, and K''{q) = 1, which yields the definition H17|l . 



4 



ring > 



(^) = y^'K'^) 



dk 

(2^ 



no) , j_± 
2/3 ^ 2/?az/ 



dk 



1 - 



1 



(27r) 



log[l + l>(fc)p[,(z.)] 



(20) 



The dominant correction to the mean field form po of the function F is 
represented in equation (jHl) by a term of the order 7'^ involving the sum of 
the tree and of the ring contributions. Adding up (fTIHl and pn|l we eventually 
find that the function F^^\u) takes the form of the derivative 



.(1) 



mg 



du 



where 



V{0) 



dk 



:2i) 



(22) 



Having derived the form of F^^^u) we can analyze now the thermodynamic 
properties of the Bose gas including the corrections of the order 7^. 



3 Equation of state beyond mean-field 
3.1 Pressure at order 7^ 

The grand canonical pressure -P(/i) satisfies the thermodynamic relation 



dfJL 



p(/i) (23) 



Considering the chemical potential p as a function of density p = /i(p) we 
get from (PHjl the formula 

The 7-expansion of the chemical potential reads 

p(p) = iimf{p) + l'^^^'\p) + 0(7') (25) 

where 

Pmf{p) = Po(p) + ap 

and po(p) = dfo{p)/dp is the chemical potential of the perfect gas (compare 
with equation (jH))). Up to the terms of order 7'^ the pressure is thus given by 

Pip) = Po(p) + ^p' + l'P^'\p) + 0(7=^) (26) 
5 



where 

p..(,)=f,,,MM ,2-) 







The correction /i*^^^(p) can be readily determined from the self-consistent 
equation Indeed, up to terms of order 7^ 

P = Po [Pm/ (p) -ap + 7'p^'^ (p)] + (Pm/ (P) - «P) (28) 

As Pmfip) — ap = po(p), upon further expanding of the first term on the 
right-hand side of (j2Hl) one finds 

p = po[/io(p)] + 7'[Po(Po(p)) P^'Hp) + i^'(Po(p))] (29) 
The identity po(po(p)) = P implies thus the relation 

F(^)(po(p)) 
Po(Po(p)) 

Moreover, since F^^\p) is the derivative of the function g(p) defined in 
we have 



^Po(p) 



dgjp) 
dp 



J M=/^o{p) 



dp 



dp 

Upon inserting (jHlj) into ()27j) we get the formula 

= .^!«) = + 1) ,32) 

where the second equality follows from integration by parts and the fact that 
5'(Po(p))|p=o = 9{p)\^l=-oo = 0. Using the explicit form (|^ of g{p) we arrive 
at the final formula 



+ V^(A;)xo(p)] (33) 



where ^ 

Xo(p) = Po(Po(p)) = P (y^^ 

denotes the compressibility of the perfect Bose gas. 

Equation ()33j) involves not only the thermodynamic functions character- 
izing the perfect gas but also the shape of the binary potential. It permits 
thus to write down the equation of state (PBj) of the interacting Bose gas 
beyond the mean-field theory. Remarkably enough, the additional pressure 
term P^^^ (p) representing the effect of fluctuations around the mean field has 
exactly the same structure as the analogous term derived by Hemmer for 
classical fluids (see H eq. (55)). 
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3.2 Low density behaviour 



It is interesting to analyze in more detail the formula (j33|l in the low density 
limit, where the compressibility Xo{p) approaches zero 



Xo(p)=/3p+f(vrATV + -, 







(34) 



The expansion of the logarithm in ()33j) yields 



Inserting here (jMj) we find 



PW(P) 



Pi 
4 



V{0)xo{p) - 2 



dk 

(2^ 



■[v{k)xo{p)r+ 



(35) 



dk 

(2^ 



[V{k)f 



(36) 



The term proportional to ~ reflects the effect of quantum statis- 

tics whereas the term proportional to (3 = l/ksT is of classical type, not 
involvoing the Planck constant. The repulsive potential V{0) > tends to 
increase the pressure, but the negative classical term acts in the opposite 
direction. Clearly, the lower the temperature the more important is the Bose 
statistics. 

Notice that the low density equation of state maintains the mean-field 
form 

Pip) = Po(p) + lay (37) 



but now with a 7-dependent constant 

1 



a + 7 



V(0)vr3/2A3 - p 



dk 



Finally, equation (jHTjl should be supplemented with the low density expansion 
of the perfect gas pressure 



Poip) 



(3^ 



Then equation (jHTj) yields the second virial coefficient at the order 7^ 
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3.3 Critical region 



At this point we make the working hypothesis that the vahdity of the formula 
(jnSl) for P^^^ (p) extends from low density up to the critical density po,c of the 
free gas provided that 7 is small enough when p is close to po,c- This is a 
plausible assumption if in the range < p < po,c, at fixed temperature T, the 
system does not undergo other phase transitions (solidification, liquefaction). 
At po,c, the compressibility of the free gas is known to diverge as 

Xo(p) ^ 



where c 



1.086 /3p2^, [5. 
p(i)(p) ^ 



Po,c - P 

For p close to po,c we can write 



(39) 




mpo,c - p) 

1 r dk 

23 J ^3 



log 



(2vr)3po,c-p + cF(A;) 
cV{k) 



1 + 



P0,c - P, 



(40) 



(41) 



2f]J (27r)= 

One sees that P^^\p) — 00 as p ^ po,c, the positive term pn|) being the 
most divergent one. Thus, choosing p close to po,c and 7^ sufficiently small, 
the pressure correction 'j^P^^\p) to the mean field can be made positive. 
We conclude that fiuctuations beyond mean field always tend to increase the 
pressure in the vicinity of po,c- 

The nature of the divergence depends on the behaviour of V{k) as /c — > 00. 
As an example we consider a power-like decay V{k) ~ bk~^, A; — >• 00, 77 > 3. 
Then the strongest divergence comes from the large values of k in the integral 
(jinj. Choosing ko sufficiently large we find 



Po,, 



2(3{po,c - p) -Jko 



~ dk 



be 



1 



27r2 
\ Po, 



k''ipo,c 



V(P0,c-p)l+3/'' 

Since the second integral (pT|l behaves as (po,c 
pressure correction diverges as (po,c — p)^(-'^+'^/^) 
corresponding to the choice 

V{r) = -^^exp{-r/ro), V{k) = ■ 



duu 



p) + bc 

be 



(42) 



u^i + be 
p)-(3/»?)^ we see that the 



Po,c 



7 7r(A/ro) 
is presented in Fig.l. 



(1 + Pr2)2 
a/3/(27rA2)3/2 = Q.l 



An illustration 



(43) 
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0.5 1 1.5 2 

Figure 1: Dimensionless pressure 1)26^ : p = P(27rA^)^/^/3 as function of den- 
sity n = p(27rA^)^/^ for the binary interaction ()43|1 . The pressure correction 
diverges at n = no,c = 2.612. The dashed hne represents the mean-field 
prediction which underestimates the pressure. 

4 One-body density matrix at order 7^ 

The grand canonical one-body density matrix R for a Bose gas in volume A 
is defined by its configurational matrix elements 

(r|i?(/i)|r') = 

'^'A 7v=l — ij ! -/A JA 

(44) 

with -f/^AT^A the A^-particle Hamiltonian and the matrix elements taken in the 
space of symmetrized wave functions. In the low density regime it admits 
the classical-like representation (equations and (ffT|l of the Appendix) 
in terms of an open polymer immersed in a gas of closed polymers. In 
the thermodynamic limit Mayer graphs resummations enable to determine 
(r|i?(/i)|r') from the formula (see Appendix, eq. (fTSj) ) 

(r|i?(/i)|r') = Fr,rK/i-ap(/i)) (45) 

where, as in 1^, Fj. j.'{i^) is the sum of multiply connected graphs with a 
root point labelled by an open polymer. It reduces to F[i>) as r = r', so 
the diagonal part (r|i?(yu)|r) = p{fi) satisfies the self-consistent equation 
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Because of the invariance under translations and under rotations Fj-ylu) = 
F(|r — r'l, I/), so that (r|i?(/i)|r') = i?(|r — r'|, /x) depend only on the distance 
|r — r'|. Hence we can set r' = without loss of generality. The Fourier 
transform 

R{p,lj,) = J dre'^ ''R{r,iJ,), p = \p\ (46) 

gives the distribution of momentum hp of an interacting Bose particle. It is 
determined by equation (j45|l in the Fourier representation 



Rip,fj,) = Fip,ix-ap{ix)) (47) 
4.1 Grand-canonical density 

We first consider the diagonal part of the density matrix (i.e. the particle 
density) up to terms of order 7'^ 

p(/i) = p„/(/i) + 7V'n/^) (48) 

where pmf is the solution of equation (jH)). Using Q and the 7-expansion (jH} 
of function F'^ we thus find 

Pm/(p) + 7V'H/^) = Po[/i-apm/(/i)-a7V'H/^)]+7'^^'^[/^-apm/(/i)] (49) 

= Po[p - apmfifi)] - i^ap^^\p)pQ[p - apmfip)] + j^F'^^^lp - ap„/(p)] 
Owing to equation (jH} the above relation yields the formula 

p^^)(p) = /^'[f-"^'"^^f, (50) 

where F^^^ is given by ()2H) . 

The density correction p^^^ (p) involves functions of the argument 

u{p) = p~ apmfip) (51) 

Notice that on one hand 

and on the other hand the mean field equation (jH)) implies 

dp p'oiu) 
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so that 

_ 1 

^ ~ 1 + ap'^{u) ^^^^ 
The density correction ()50|) can be thus rewritten in the form 

p«(/.) = F«[K^)]|^ (53) 
which in view of the structure of equation (PT|) finally yields 



where function g{iy) has been defined in 

This correction to the grand-canonical density can also be studied at 
low density and in the critical region. If z/ ~ /i — — oo (i.e. low density) 
one finds that p^^^i^) ~ A{P)e^l^^ where, as for the pressure, the sign of 
the coefficient A{(3) depends if effects of Bose statistics dominate classical 
corrections or not. Extrapolating the formula ()54|) to the neighborhood of 
the critical mean field chemical potential fimf,c = CLPo,c, one deduces first 
from the mean field equation that z/(/i) ~ —C{P){fimf,c — fJ'Y,C{P) > 0, 
as /i — > fimf,c{fJ' < l^mf,c)- Then analyzing the k-integrals as in subsection 
3.3 one sees that p^^^p) diverges to — oo as /i — Prnf,c (for V{k) ~ 
p^^\p) ~ —{Pmf,c — p)~^^~^^^'^'^^ as for the pressure). Thus, taking p close to 
Pmf,c and 7 small enough, the density decreases when fiuctuations are taken 
into account. This is illustrated in Fig. 2. 



4.2 The momentum distribution 



The momentum distribution R{p, p) at order 7'^, considered as a function of 
the density 

R{p, p) = m (p, p) + (p, p) (55) 

is found from the relation (jTTj) by using the 7-expansion of F"* which follows 
from dHSl) an from the relation u{prnf{p)) = [Pmf - «p(/^)]m=m,„/(p) = Po(p)- 
We observe that at the dominant order (see (f7njl . (f77|l ) 



(56) 



is the free Bose momentum distribution at density p. Proceeding as in ()28|) - 
(jHUjl one finds 



M'\p, p) = F('\p, poip)) - 



n 



Po(^j 



(57) 
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Figure 2: Dimensionless density ()49|) : n = p(27rA^)^/^ as function of the 
chemical potential m = for the binary interaction (j43|l . The dashed line 
represents the mean-field prediction which overestimates the density. We put 
here a = a/3/(27rA^)^/^ = 1.5, a(7A/ro)^ = 4y^2/7r. The density correction 
diverges at m = mo,c = opo,c = 3.93. 



Inserting the expression ()85j) for F^^^'^p, v) and for F*^^)(z/) = / dpF^^\p, v) / (27r) 
one obtains the final result 

where 



J v=no{p) 



c(p,/?) 



2/?y (27r)3 Vl + n%o(/^o(p)), 



Notice that 



as requested by the fact that the integrals 



(ip 

(2^ 



(5^ 



(59) 



(60) 



are fixed by the total density. 

Some properties of -R(p, p) can be derived from those of the free Bose 
distribution ^ 

^o(p, z) = 

exp(A^p^/2) — z 
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where z = e^'^, < 2; < 1, is the activity parameter. One sets 

ni{p, z) = p'^n^ip, u) = no{p, z){l + no{p, z)) 

n2ip,z) = P~'^n'^{p,u) = ni{p,z){l + 2no{p,z)) (61) 



Then 
with 



R{p, z) = no{p, z) + e[ni(p, z)r{z) - n2(p, z)] (62) 

^ ^ (3Mu) Jdpn,{p,z) ^ ' 

The parameter e = 7^c(p, /?)/3^ in (jU^ incorporates all the p-independant 
factors (see eq. (j58|l ). Since e is proportional to 7^ it can be chosen as small 
as one wishes at any density p < po,c- We observe that 

n2{p, z) < ni{p, z) [l + Y^r^) = ^) (137) ^^^^ 

implying 



r{z) < ^ 

z 

[ni{j),z)r{z) - n2{p,z)] 



P=0 (I _ z)2 



r{z) 



< (65) 



One concludes that the distribution at p = in presence of the interaction is 
always less than the free value z/{l — z). Evaluating the correction asp ^ cxo 
gives 

[ni{p, z)r{z) — n2{p, z)] ~ zexp(— A^p^/2)(r(z) — 1) > (66) 

since no{p,z) ~ ni{p,z) ~ n2{p,z) ~ zexp(— A^p^/2), p ^ 00, and for all p 
and z, n2{p,z) > ni{p,z), implying r{z) > 1. Thus the momentum distri- 
bution is flattened and broadened by repulsive interactions, as illustrated in 
Fig.3. 



5 Concluding remarks 

One can make the following points: 

(i) Divergences occuring at the critical values of the free gas reflect the 
fact that the order 7^ of the correction is not adequate there. Since we know 
from rigorous work that the exact pressure converges to the mean field 
pressure Pmf for all values of the thermodynamical parameters jB] one may 
conclude that the rate of approach to mean field at critical values of the 
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free gas is of the order 7'' for some < 5 < 3. The determination of this 
exponent is an open problem which necessitates a more thorough analysis of 
Mayer series. 

(ii) In our analysis, the critical values of density and of temperature are 
still those of the free gas. One could expect that for non-zero but small 
7 the critical mean field chemical potential /im/,c = «Po,c and density po,c 
will be sligthly displaced. Such an information cannot been obtained here, 
but requires, for instance, a deeper non perturbative understanding of the 
behaviour of the partition of a single polymer |TT| as also discussed in Section 
V.A of p. 

(iii) The results of the paper enable nevertheless to explore the imme- 
diate vicinity of the critical mean field values by taking 7 sufficiently small 
(provided that our formulae keep their validity up to this point). Then qual- 
itative statement could be established, as the fact that the exact pressure 
should be higher than that predicted in this region by the mean field . An- 
other conclusion from our analysis is that approaching the critical mean field 
chemical potential /im/,c = apo,c along an isotherm one finds the interacting 
gas at a density lower than that appearing in the mean field approach. 
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A Appendix 

A.l Polymer representation of the one-body reduced 
density matrix 

The representation of the off-diagonal reduced density matrix (r|i?(/i)|r') in 
the formalism of random polymers described in Section 2 of PP requires the 
consideration of open Brownian paths (open polymers) in the time interval 
< s < q with extremities at r and r'. In the Brownian bridge notation they 
are parametrized as 

1 - - 1 r+ -r' + AX(s), < s < g. (67) 
qj q 

where X(s) is a closed path distributed with the Brownian bridge measure 
(eq. (18) in [T] ; here and in the sequel we use the same notations as in 
[l]). The open polymer £r,r' = (r, i"', <?, X(s)), < s < g is characterized 
by its end points r, r', the number q of particles belonging to it and its 
random shape X(s). By a slight generalisation of the analysis that led to 
the "magic formula" (eq. (14) in ^J) one obtains the density for an open 
polymer immersed in a grand canonical ensemble of closed polymers as 

Pop(-^r,r') 

~ E 7 TTT^('^r,r') / IT dCiZ{Ci) exp(-^f/(£r,r' , -^2, • • • , -^n))- 

(68) 

where Sa is the partition function. If r = r' this reduces to the the loop 
density pioop studied in ^. In view of its form (jUHjl Pop('^r,r') has the 
Mayer expansion presented in eqs.(29)-(31) of JLJ in terms of Ursell functions. 
The only difference is that the argument of the root point 1 = £r,r' has to 
be identified with the open polymer (jHTj) . Polymer interactions, effective 
activities and Mayer bonds are the same as in jl] except at the root point, 
where the interaction between a loop Cj = (Rj, qj, Xj) with an open polymer 

^The quantity — ksT log pop^C^.r') can be interpreted as the excess grand potential 
when a fixed open polymer £r,r' is introduced in the system of closed polymers. 
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IS 



V{Cry, Cj) — 

£ ds £ dsjS{s - Sj)V N 1 - - j r + -r' + AX(s) - R^- - XXj{s 



(69) 

and the self-energy of the open polymer is 



U{Cr,r') = ^J^'dsiJ^'ds2S{si-S2)v(^^^-^{r-r') + X{X{s,)-X{s2))j 
- IdViO) (70) 

Finally, to obtain the one-body density matrix {r\R{^)\r') , one has to inte- 
grate on the internal degrees of freedom of the open polymer 

(r|i?(/i)|r') = ^qj D,(X)exp (--^-^) Po^^ry) (71) 



q=l 



The q factor takes into account the presence of q particles in the open polymer 
and the additionnal Gaussian comes from the Brownian Wiener weight for a 
path starting in r at time s = and ending at r' at time s = q. From now 
on the analysis of the Mayer series can be performed exactly along the same 
lines as that given in ^ with the following results. 

(i) The Mayer series representing (r|i?(/x)|r') converges for /i sufficiently neg- 
ative (yU < — apo,c); for a scaled potential V^(r), the convergence is uniform 
with respect to 7. 

(ii) Let /(£r,r') be the value of the sum of all multiply connected graphs with 
root point Cj.y, and define the function of the chemical potential u 

FrA^) = Y.qJ D,(X)exp (^-^^J ^^r.r')- (72) 

Then the reduced density matrix is given by 

(r|i?(/x)|r') = Fr,rK/i-ap(/i)). (73) 
where the density p{fi) solves the self-consistent equation (jS)). 



16 



A. 2 The 7^ correction 



For a scaled potential (7 <^ 1), the function F''{r,h') takes the asymptotic 
form 

F^r, u) = F°_(r, u) + 7'[F,L(r, u) + F,\„^(r, u)] + 0(7^) (74) 

as in (jni) where F/°ee(r, z/) represents the sum of zero order contributions 
coming from the tree diagrams. The sum of tree diagrams yields 



,3/2 



K,'''{q), quadr 



(75) 



so that 



.(0) 



tree V ' 



r, z/ 



lim F7 

7^0 



1 °° exp 

tree[r, v) = (2vrA2)3/2 



2A2q 



.3/2 



(2^ 



-ikr 



no{k, v) 



(76) 



is nothing else than the off-diagonal reduced density matrix po(r, = (r|i?o('^) |0) 
of the free gas at chemical potential v . Indeed, F/°gg(r, z/) is the Fourier 
transform of the Bose occupation number density no{k, u) 



no{k, v) 



1 



exp 



(A|fc|)^ 



-(3v 



(77) 



The 7^ correction is obtained in expanding the partition function hC[q) of a 
single polymer in FjJ^^^^r, z/) as done in eqs (fTT|) to (fTBjl 



tree V ) 



r, u 



WO) 
2 

1 d'' 



(/«(r,z/)-/(2)(r,z/) 
Po(r, J^)- 



with 



(7J 



Clearly f^''\T,u)\r=o = f^^\v) defined in We now perform the ring 



summation to determine Fj:llg(r,i') . The contribution to Iringi^ry) of a 
ring with one root point £r,r', n integrated vertices Cj and n + 1 linarized 
bonds {—PV) is 

l[z{q){-PV{q,Cj+l)) hPV{Cr.,Cr,r')) (79) 
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The factor 1/2 is the symmetry factor of the graph. Introducing the scaled 
potential j^Vl'jr) and changing the n spatial integration variables Rj to 7Rj 
will produce an overall factor 7^ in ()79|1 since there are n+1 bonds. At order 
7'^, we are thus entitled to neglect the quantum fluctuation part 7AX(sj) 
in the arguments of the potential (use dominated convergence and the fact 
that the Gaussian measures Dq.i^j) are normalized). For the same reason 
we neglect the 7 dependence in the activities. Within this approximation we 
take the bonds and the vertices equal to 



q{2Ti\'^qfl^ 



z'-'\q) 



^0) 



At this point we keep the 7 parameter in combinations 7r and 7r' since r, r' 
can be large. Hence after the change of variables Rj — > 7Rj ()79p becomes 



2!.(o) 

2^ 



(q) JdRi-'-dRn -pjysVni- 



7r H — 7r 
1 



Ri 



s 
Q 



7r 7r 

q 



Introducing the Fourier transform V{k) and using the convolution theorem 
yields 



^;2(°)(g)(/«(z/))" fds r ds' 



dk 

(2^ 



exp 



s — s 



n+1 



where f^^\i^) = Y.'^i l"^ ^^^Hl) is the function (fTTj) . This has to be summed 
on n = 1, 2, . . . and integrated over the internal degrees of freedom of the 
root point according to (f7^ . One finds 



-^g2;(°)(g)exp 



r — r 



(27r 



i'j k 



2gA2 J 
f s-s' 
V 1 



J l+pV{k)fWiu) 



(83) 



(82) 
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Finally, taking (jHn|l into account, the Fourier transform F^^\p, ly) = J dre'^^''^Fj:llg{r, v) 
at wave number p is found to be 



exp 



A2 / 'X 2 



2 



. I s - s 
p + ^yk 



(27r)3 



l + PV{k)f(^){u) 



As a consequence of the Gaussian weight occuring for off-diagonal elements 
r 7^ r' the remaining 7 dependance has been neglected in fl84|) . p'^\p,i^) is 
the Fourier transform of the function (f75|) . We note that 



/ 



(2vr) 



3 ring\ri " j rmg 



reduces to the formula (fTHj) used in Section 2. 

The final result for the 7^ correction is obtained by adding the tree con- 
tribution dZHl) to dHl) 

F^'\p,u) = FBL{p,i^)+Fillg{p,u) 

- (35) 



2 ^ 2/? 'J (2vr)3 l + V(A;)p^(z. 

where we have expressed all quantities in terms of the density and the Bose 
momentum distribution ()77|1 of the free gas. 
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